Influence of initial size on higher cumulant ratios of net-proton number fluctuations 
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With the help of AMPT default model, we study the influence of initial size (centrality of collisions) 
on higher cumulant ratios of net-proton distributions. If the centrality is presented by impact 
parameter, there is a strong centrality dependent, in particular, in those peripheral collisions. This 
dependence is slightly reduced if the centrality is presented by number of participant, or charged 
multiplicity. However, the dynamical ratios are almost centrality independent. So the centrality 
dependence of dynamical ratios at RHIC beam energy scan are presented. 

PACS numbers: 25.75.Nq, 12.38.Mh, 21.65. Qr 



I. INTRODUCTION 

One of the main goals of relativistic heavy ion collisions 
is to map the QCD phase diagram in particular, the 
critical point. Recently, the higher cumulants ratios of 
net-proton distribution are suggested to be a good probe 
of QCD critical point (QCP) [2|-|6||. They are skewness 
and kurtosis defined as, 



S = 



((57V 3 ) 
(SN 2 ) 3 / 2 ' 



(oY 4 ) 
(SN 2 ) 2 



-3. 



(1) 



Where SN = N — (N), and N is the number of net- 
proton. The skewness and the kurtosis present respec- 
tively the asymmetry and the peak of the distributions. 
For a Gaussian distribution, they are zero. While near 
the QCD critical point, a non-Gaussian distribution is 
expected. They deviate from zero, and change the 
sign 0,i]. 

However, it is found afterwords that those ratios are 
influenced by various trivial effects, such as the central- 
ity bin width Q, and statistical fluctuations due to the 
finite number of proton produced at RHIC incident ener- 
gies [1, Q . To eliminate influence of centrality bin width, 
the centrality weighted method are recommended ■ To 
subtract the statistical fluctuations, the dynamic ratios 
are suggested, 
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Where statistical parts of proton and anti-proton are 
determined by Poisson-like fluctuations. The statistical 
fluctuations of net-proton is given by a Skellam-like dis- 
tribution. So, S s tat = and « stat = -pjy. Here, M 
is the mean of total proton number. They are both ze- 
ros if proton and anti-proton are emitted as independent 
Poisson. 

In addition to the above mentioned effects, the initial 
size, i.e., the centrality of collisions, should also influence 
the measurements. The centrality of collisions is origi- 
nally defined by impact parameter (6). It is the distance 
of two incident nuclei. It takes minimum and maximum 
in central and most peripheral collisions, respectively 
However, it is not experimentally measurable. In experi- 
ment, the number of participating nucleons (Y part ), and 



charged multiplicity (Y c h) [HI are usually used to quan- 
tify the centrality. In fact, only charged multiplicity of 
final state is directly measurable. From any one of these 
three variables, the other two variables can be obtained 
by the Monte Carlo simulation of Glaubc model. Conven- 
tionally, the centrality in model calculations are classified 
according to the impact parameter. While in real data 
analysis, it is usually given by the number of participants, 
which is deduced from directly measured multiplicity. 

The influence of centrality definitions varies from ob- 
servable to observable. Some observables are not sensi- 
tive to the definitions, like elliptic flow V2- But others are 
not, such as forward-backward correlations, and system 
size. Their centrality dependences change significantly 
if the definition of centrality changes from one to an- 
other 1T|. In this paper, we will study the influence of 



centrality definitions on higher cumulant ratios using the 
AMPT default model 0. 



II. INFLUENCE OF CENTRALITY 
DEFINITIONS 

In AMPT default model [lH, we can easily classify 
the centrality by impact parameter, number of partici- 
pants, and multiplicity, respectively. They are usually 
quantified by percentile. So in Fig. 1, we present the 
centrality dependence of Kurtosis (a), dynamical kur- 
tosis (b), sa(c), and k<j 2 (d) for Au + Au collisions 
at 39GeV, where ccntralitics arc defined by impact pa- 
rameter (red open stars), number of participants (blue 
open crosses), and multiplicity (black open triangles) in 
\r/\ < 0.5, respectively. All ratios are calculated in the 
same phase space cuts as RHIC/STAR data ||. The 
centrality weighted method is also used here to eliminate 
the effects of centrality bin width. 

We can see from 4 sub-figures that red open stars 
are obviously higher than both blue open crosses and 
black open triangles, in particular, in peripheral colli- 
sions. While later two points are close to each other. 
This shows if the size fluctuations are presented by im- 
pact parameter, it contributes more to the ratios. The 
ratios arc sensitive to the initial size fluctuations. So 
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FIG. 1: Centrality dependence of kurtosis (a), dynamic kurtosis (b), S*ct (c) and K*a 2 (d) for Au + Au collisions at 39 GeV, 
where centrality are classified by Apart (blue open crosses), A cri (black open triangles), and b(red open stars), respectively. 
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FIG. 2: Centrality dependence of dynamical skewness (a) and 
kurtosis (b) for Au + Au collisions at 7 RHIC/BES energies 



We can see from Fig. 2(a) that at 2 lowest incident 
energies, the dynamical skewness are negative in periph- 
eral collisions and becomes smaller and smaller towards 
the peripheral collisions. At 5 higher incident energies, 
they are around zero and almost energy and centrality 
independent. The dynamical kurtosis in Fig. 2(b) are in- 
cident energy independent. They are all negative in pe- 
ripheral collisions, become smaller and smaller towards 
peripheral collisions, and are close to zero in mid-central 
to central collisions. So in peripheral collisions, there are 
still contributions in additional to Poisson-like statisti- 
cal fluctuations in the measurements. The other effects, 
such as cuts of phase space and so on, may contribute to 
it too, and related studies are on going. 



whenever we compare the data with model calculations, 
the same centrality definition should be used. 

It is interesting to notice that centrality dependence of 
dynamical kurtosis in Fig. 1(b) are well reduced if the 
centrality is defined by number of participants, or mul- 
tiplicity. Moreover, they are close to zero in all central- 
ities, in contrary to the kurtosis given in Fig 1(a). So 
dynamical ratios well reduce the influence of the initial 
size fluctuations. 



III. ENERGY AND CENTRALITY 
DEPENDENCES OF DYNAMICAL RATIOS 

As we know, there is no critical fluctuations in AMPT 
default model. The results from the model can serve as 
a good reference of background effects. So the energy 
and centrality dependences of dynamical skewness and 
kurtosis for Au + Au collisions at seven RHIC/BES in- 
cident energies, 7.7, 11.5, 19.6, 27, 39, 62.4, 200GcV, are 
presented in Fig. 2(a) and (b), respectively. Where the 
centrality is defined by A^h, and plotted by Af part , same 
as experimental data. 



IV. SUMMARY 

With the help of the AMPT default model, we study 
how the initial size (centrality of collisions) influence the 
higher cumulant ratios of net-proton. It is found that if 
the centrality are presented by number of participant and 
multiplicity, their corresponding centrality dependence 
are similar. But they are quite different from that which 
centrality is presented by impact parameter. So when we 
compare the experimental data with model calculations, 
the same centrality definition should be used. It is also 
found that dynamical ratios well reduce the influence of 
initial size. So the centrality dependence of dynamical 
skewness and kurtosis for Au + Au collisions at seven 
RHIC/BES energies are presented. They can served as a 
good reference of background effects. 
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